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Continual interest in quest for perfeet imaging in optieal design, reeently reinvigorates the 
efforts to advanee the imaging resolution to eater for human’s visual observation or other 
entertainment requirements. Despite the material-loss-suffering eoncept of perfect lensing 
via negative refractive index [1, 2, 3, 4, 5], an alternative approach to achieve subwavelength 
imaging using positive refractive index was proposed [6], namely mirrored Maxwell fisheye(a 
finite spatial-varying profile), which was based on manipulating curved geometry of light [7]. 
This idea to design the geometry for light de facto gave birth to a novel sub-field of 
transformation optics [8, 9, 10, 11, 12]. A worth-noticing paper [13] declares that an active 
drain in a mirrored Maxwell fisheye does not give subwavelength resolution of imaging. 
Consequently, an extension to combine coordinate transformation(herein precisely conformal 
mapping) and spatial profile(herein mirrored Maxwell fisheye profile) showed an interesting 
duality of simultaneous cloaking and imaging [14], but no thorough investigation on the 
imaging resolution was performed. This manuscript’s novelty is two fold: one in treating 
this resolution aspect in a transformed full Maxwell fisheye, the other in shrinking an infinite 
Maxwell fisheye profile into a finite unit circle by conformal map [8, 12]. The previous 
investigation [14] made use of a mirrored Maxwell fisheye, which necessitated a forbidden 
cavity region delineated by a curved shell of perfect electric conductor(PEC). An equivalent 
profile therefore, can also be investigated: the, full Maxwell fisheye without this unnecessary 
shell of PEC. This equivalence relies on the fact that their optical path between the full and 
the mirrored Maxwell fisheye are equal [6] . Here, we propose another heterogeneous medium 
for light to perform similar functionality of good imaging, as proposed in [14]. 


2. Profile design for proposed imaging 


Eet us start from the full Maxwell fisheye profile for refractive index. 


n{r) = 


l+(f)^ 


,r 6 [0, oo) 


( 1 ) 


where all light rays make circles(rn, is a characteristic radius). It was Euneberg who realised 
that Maxwell fisheye profile can be generated by stereographic projection from a sphere 
surface to a plane [10, 11]. By the idea that light travels along geodesics along sphere(in 
virtual space), an optical implementation of imaging device follows: all rays emitted from 
any point in full Maxwell fisheye will converge at its conjugate point-antipole-to form an 
image point [6](cf. converging circles in figure 1(a)). An obvious fault for the full profile 
of 2D Maxwell fisheye is its infiniteness in area and hence the truncation trick along with 
mirror was employed [6]. However, here we attempt to sidestep this obstacle by applying 
conformal-map method [8]. 

We use the conformal-map to transform imaging behaviour of the full Maxwell fisheye. 
In principle, one can implement any conformal map to this Maxwell fisheye profile to create 
a new heterogeneous medium to perform similar imaging functionality as in figure 1. A 
conformal map creates a virtual space where light travels but this virtual space is made of a 
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stack of Riemann sheets in general [8, 15]. In this paper, Zhukovski map is chosen which 
creates two Riemann sheets which are connected by their branch cut [8, 14]. Because the 
full and mirrored Maxwell fisheye preserve the same optical path(cf. virtual sphere in Figure 
3(A) [6]), the mirror in previous work [14] can be removed so as to dismantle the constraint 
of the curved shell of PEC, by simply taking the full fisheye on the lower sheet of virtual 
space. Equivalent to ray-tracing in figure 1(a), one obtains a wave field picture Efw) [6] on 
the lower Riemann sheet (half of the virtual uniform space) in figure 1(b). 


2.1. Conformally-mapped Maxwell fisheye 


For a fair comparison with our previously-published wave picture [14], we stick to Zhukovski 
map 

cfi 

w{z) = z -I- —, (a > 0) (2) 

z 

in optical conformal mapping method [8] although principally any conformal map reserves 
the imaging behaviour of Maxwell fisheye (1). Under Zhukovski map, a physical space(plane) 
is mapped onto two Riemann sheets, of which the lower and upper sheets are geometrically 
connected by the branch cut-a segment from -2a to 2a. The two sheets are labelled by two 
Cartesian coordinates {u, v) and {u', V) shown in figure 2(a). 

We choose to shift the centre of Maxwell fisheye profile, M^, away from the origin by an 
arbitrary distance b on lower sheet. The merit of using conformal map in this case is to transfer 
imaging point from underneath to above the ground. Therefore one has to put the image 
point(marked in hollow circle on the lower sheet in figure 2(a)) on the branch cut [-2a, 2a] 
(illustrated by purple segment in figure 2(a-b)) and thus the source(marked in blue cross in 
figure 2(a)) away from that. This is because a pair of image and source in a full fisheye appear 
on opposite sides of the centre(marked in red dot) along a line joining the source and the 
image(cf. figures 1 and 2(a)). If the source point fell on the branch cut, its equivalent in z 
physical space would degenerate into two sources(which images to two points separately) off 
w'-axis on the unit circle(branch cut), when parameter b is swept over under the condition 
rw < a. In this case, the image point would instead lie underground(in the unit circle) in 
the virtual space and one feature of conformal imaging [16]-the double degeneracy of one 
source point imaging to two points- is broken. Furthermore, the ray-tracing computation in 
forthcoming section 4 is doomed as well. This falls beyond the scope of this manuscript and 
we restrict the fisheye size under the condition r^ > a. 

Thus based on the optical conformal mapping method [8, 11], the refractive index 
profile(in physical space, its range ~ [0, 33]§) for our conformally-mapped full Maxwell 
fisheye is written as 


n{x,y) = n 


dw(z) 


dw(z) 

dz 


dz 


1 , 


2 

w-b |2 _|_ ’ 


for |z| > a; 


for |z| < a. 


z = x-\-iy 


( 3 ) 


§ It is interesting to notice that the zero points of refractive index occur at z = +1, which aligns with the doom 
of trajectory-tracing from Hamilton’s equations at these two points. 
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Figure 1. (Color online) (a) Ray-tracing diagram for a full Maxwell fisheye(focus (0.5,0) 
originated from source (-2,0)), in which red dashed circle delineate at r = in a full 
Maxwell fisheye; (b) Wave value distribution Edw) in arbitrary unit corresponding to (a) 
according to [7]. Parameters for (a-b); image point(0.5,0), scaling factor a - 1,fisheye 
radiusrvv - \-,l - 2A,y - arctan2,;^f - n in (3), (4) in the forthcoming section 2.1. 


Considering TE polarization where eleetrie field direets perpendieular to the eomplex z plane, 
the sealar Helmoltz equation stands for any wavelength [6]: 


A = 


2nr„ 

sjl{l + 1 ) 


,/> 0 . 


( 4 ) 
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To compare fairly with previous simulation [14], we follow to take integer values of I in this 
manuseript but I ean be non-integer as well, whieh would vary the imaging phase. 

Notice that the full Maxwell fisheye profile bends the uniform eoordinates (m',v' 
Cartesian)on lower w sheet into curvilinear ones (0, (p spherical) as plotted in lightblue and 
lightgreen in figure 2(a), according to inverse stereographie projeetion(One can calculate arc 
tangent value via the eommand of ArcTan[u’-b, v’] in Mathematica 10.0 [17], taking 
into aeeount with whieh quadrant the point (u' - b, v') is in). 



( 5 ) 


(p = arotan(M' - b, V) 6 [0, lit). 


( 6 ) 


An imaging proeess eonsists of the light flow from any souree point(marked in blue cross) 
towards its antipodal(in hollow circle) and then leakage(refraetion) onto the upper sheet, 
illustrated by three representative light paths in thiek green eurves in figure 2(a). It is 
interesting to notiee that all rays emanating from the souree do not refraet on the upper sheet 
since total refleetion eould oeeur for partial of them. Thus the light trajeetories on the lower 
sheet only occupy a finite area, out of whieh is the forbidden zone for light. 

3. Wave simulation 

In two-dimensional ease, put an aetive point souroe(lA) at 5(z) on z plane, we simulate 
on the TE field distribution(£' 2 (x,y)) via COMSOL Multiphysics [18]. As a preliminary 
demonstration, we plot a distribution of TE field in figure 3 under a certain set of parameters. 
Notice a truncation is taken beyond certain radius away from unit circle(for instance, 10r„,) 
sinee the heterogeneous medium there ean be approximated as unity similar to [14]. We 
ehoose to position the souree so that the eorresponding image falls at the middle of the braneh 
eut-the origin in virtual spaoe(+i in physieal spaee), as seetion 2 explains. 

This feature of unity refraetive index far away from the unit cirele, allows one to position 
two image detector in the free spaee environment, which may be easy to realise in experiment. 
The diseontinuous pattern at the braneh eirele |z| = 1 ean be attributed to the slight jump 
between their refraetive indiees of two Riemann sheets. This diseontinuity is also present in 
previous simulation work [14](Eig. 2 therein), whieh may be improvable by designing a layer 
of gradient index smoothly varying but we stick to the simple profile (3) for simplicity. In 
figure 3(b), the light originates from source point S{z) inside the unit eirele and diverges out 
on two distinetive paths upward until to the image point i on the upper half plane and -i on the 
lower half plane as well(mirror-symmetry), similar to the ease of mapped mirrored Maxwell 
fisheye [14]. 

In order to seek the best resolution of imaging, we sweep over several parameters to 
investigate on the image resolution with respect to the operating wavelength in aeeordence 
with (4). 
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(b) 



Figure 2. (Colour online) (a) Illustration of the lower (m', v') and the upper (m, v) Riemann 
sheet. On the lower sheet, red dashed circle indicates the coordinate r„ - I, blue contours 
for 0 - arccosz/r^, 6 [0, tt) and green for (p 6 [0,27r) as curved coordinates. The red 
dot indicates the centre of the Maxwell fisheye put on the lower sheet. Thick green 
circles(including a straight line through the branch cut in purple) are three light trajectories 
contributing to form an image at the origin point in hollow black circle 0„ from a source point 
in blue cross S„: - b,M„S„ - t^jb. Notice that the light trajectories deflect at the 

branch cut and transmit onto the upper sheet(reflection overlooked). The wave distribution 
Eiu') [6] and the refractive index n{u') for the branch cut on the lower sheet are plotted on the 
upper sheet. Parameters: b - \,r„ - A,l - 24. (b) Virtual curved space (0,0) on a sphere, for 
the lower sheet in (a) under stereographic projection. On the sphere light travels in geodesics, 
in green circle, diverging from blue source and converging at black image point(geodesics 
equation based on Eq. (19.9) [11] with the corresponding initial conditions). 


3.1. Sweeping the centre of the fisheye via b 

Firstly, to find an optimal image resolution for the considered profile, we opt to shift the 
centre of the fisheye along M'-axis on the lower w sheet(sweeping b). Due to the geometrical 
symmetry of the case, varying b within positive regime is sufficient to disclose the resolution 
scaling. In figure 4(a), average energy flowing out along the upper half circle(branch cut) 
is plotted versus radiant angle, under different values of b. It is found that as b decreases to 
zero, the image resolution becomes sharper until a lower limit of around half wavelength if we 
define the imaging resolution by full width at half maximum(FWHM) as shown in figure 4(a) 
inset. 
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(a) Surface: Electric field, z component emw.Ez (V/m) 


(b) Surface; Electric field, z component emw.Ez (V/m) 




Figure 3. (Colour online) (a)Electric field distribution (in V/m) in z direction(COMSOL 
simulation); (b)Zoom-in of (a), (a-b) Unit circle is marked black to indicate the boundary 
between two Riemann sheets in virtual space. Parameters; b - \,r„ - A,l - 2A. 


3.2. Sweeping the size of the fisheye 

Secondly, the geometrical size of Maxwell fisheye on the lower sheet can be varied as well. 
Although the imaging resolution appears shrinking as increases, it never shrinks below 0.3 
wavelength in figure 4(b)(in physical space, the actual wavelength is Aln{z = i) = All). The 
sharper trend with increasing may imply the refractive index jump along unit circle |z| = a, 
ignored by the previous study [14]. As r„ increases, the discontinuity in the refractive index 
n{z) at the unit circle diminishes and thus better resolution is achieved. This discontinuity in 
the refractive index may also explain the scattered colour pattern in figure 3. 

2.2. Varying the wavelength via eigenvalue I 

A similar trend can be observed when one tunes integer values of I in figure 4(c), that 
the imaging resolution evolves worsen as I increases but still no shaper resolution than 
0.2 wavelength(corresponding to the least value 1 = 2) can be achieved herein, which is 
comparable to the resolution obtained in [14]. As mentioned in 2.1 about (4), only integer 
valuesof I are taken but non-integer value applies as well. 

To summarise section 3, we sweep over three parameters to observe the imaging 
resolution dependence on them. The optimal parameter according to our simulation is: /? = 
0,rw = S and 1 = 2 and the optimal resolution achievable in our scenario is 0.2 ~wavelength 
according to our definition of FWHM. 

4. Ray-tracing diagrams 

The bifurcating behaviour of light wave, inside the unit circle in figure 3 tempts one to wonder 
the propagation scenario of light there. One simple method is to investigate into the ray- 
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Figure 4. (Color online) Average energy flowing out along the upper half circle when one 
sweeps over (a) b, (b) r„ and (c) 1. Note that (b-c) main panels and all three insets are all 
normalised in ordinate and abscissa to directly compare with the operating wavelength in 
physical space:/l/n(z = i) - Ajl according to (3) and (4). Parameters: - A,l - 24(the 

same wavelength A - 0.5130m as in [14]); (b) b = 4, / - 24; (c) b - 0, - 4. 


tracing diagrams. When it comes to ray-tracing in virtual space, it is straightforward to predict 
that the light rays turn out to be circles underneath and straight lines above(cf. figure 2(a), 
if one neglects reflection at the branch cut||). However, it may be of interest to peep into 
ray tracing in physical space in transformed Maxwell fisheye. This section will address this 
concern. 

Based on the refractive index distribution (3), we use Hamilton’s equations dfT = dcojd^ 
and (9fk = -d(o/dr to compute the ray trajectories [11, 20] by Mathematica [17], in which k 
indicates the momentum, r the position, co the Hamiltonian. An instance is plotted in figure 5 
for the same parameter as figure 3 to illustrate the light flow of imaging: light rays originate 
from source point z = 0.0590, diverge at two image points z = +i and then propagate outward. 
This ray-tracing diagram aligns remarkably with the wave solutions we obtained in figure 3. 

II Notice that the discussion for null reflection from an incidence wave from the upper Riemann sheet (on 
equations (11)-(13), p4 [19]) does not apply since in our scenario the incidence wave comes from the lower 
sheet instead of from the upper in [19]. 
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Ray tracing (b=1 =4) 

y/a 



Figure 5. (Colour online) Ray trajectories computed according to Hamilton’s equations. Unit 
circle is marked black to indicate the boundary between two Riemann sheets in virtual space. 
Rays starting with different directions are coloured differently for illustration. Parameters: 


5 . Summary 

In this manuscript, we demonstrate a new seheme to design heterogeneous optieal 
medium(isotropie) from perspeetives of wave simulation and analytie approximation, by 
making use of full Maxwell fisheye as well as eonformal-mapping method [8, 12]. The 
merit of this idea lies on removal of truneation of the profile of Maxwell fisheye, and that 
it maintains the same imaging resolution as a truneated Maxwell fisheye(with mirror). A 
possible alternative to eompressing infinity into finiteness as eonformal map does, eould 
be the inverse tangent transform in Penrose diagrams(ef. pp201-3 [11]). There is more 
versatile funetionality to investigate on, if one puts on the lower sheet more than one light 
souree, for instanee the mutual interferenee between sourees, whether eonstruetively or 
destruetively [14]. 
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